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Cosmological consequences of a coupling between massive neutrinos and dark energy are inves¬ 
tigated. In such models, the neutrino mass is a function of a scalar field, which plays the role of 
dark energy. The evolution of the background and cosmological perturbations are discussed. We 
find that mass-varying neutrinos can leave a significant imprint on the anisotropies in the CMB and 
even lead to a reduction of power on large angular scales. 


The discovery of the accelerated expansion of the uni¬ 
verse is a major challenge for particle physics (see 0 
for latest results). According to General Relativity, the 
dynamics of the universe is dominated by a new (dark) 
energy form with negative pressure. A well-motivated 
candidate for dark energy is a light scalar field Hi 
From the particle physics point of view, however, such 
light scalar helds are problematic: first of all, why is the 
mass so small and how can this low mass be stabilized 
against radiative corrections 0]? Secondly: is this field 
coupled to any other matter form? And if not, why not 
m? Indeed, early papers on quintessence discussed this 
possibility in detail with models in which dark matter is 
coupled to dark energy The discovery of “dark 

energy” clearly requires new physics for its explanation. 

In this work we explore the cosmological consequences 
of an idea recently put forward in [^. According to this 
idea, dark energy and neutrinos are coupled such that 
the mass of the neutrinos is a function of the scalar field 
which drives the late time accelerated expansion of the 
universe. In general, the held will evolve with time and, 
hence, the mass of the neutrinos is not constant (mass- 
varying neutrinos). One of the motivations for such con¬ 
siderations is the question of whether there is a relation 
between the neutrino mass scale and the dark energy 
scale, that has a similar order of magnitude compared 
to the detected neutrino mass splittings. In such mod¬ 
els the origin of the neutrino mass and dark energy are 
interlinked. Astrophysical and cosmological implications 
of such models have recently been studied in Q. Here 
we study for the hrst time the transition of coupled neu¬ 
trinos from the relativistic to the non-relativistic regime 
as well as the dynamics of the dark energy field. We also 
consider how the coupling affects the cosmic microwave 
background radiation (CMB) and large scale structures 
(LSS). 

The dark energy sector is described by a scalar field 
with potential energy This potential has to be 

seen as an effective, classical one, since the coupling be¬ 
tween the scalar field and the neutrinos can lead to sig- 
nihcant quantum corrections 0, a problem also present 


in models with dark matter/dark energy interaction j3|. 
To be specific, in this letter we will choose a standard 
quintessential potential, namely the exponential poten¬ 
tial = Vb exp(—-s/SA/Zv^) ( in the following we set 
8ttG = 1) 13. With this choice, our theory differs from 
the one proposed in |^. There, the choice of potential 
was such that the mass of the scalar field is much larger 
than the Hubble parameter H from times before big bang 
nucleosynthesis until today. In contrast, with our choice 
of potential, the mass of the field will be at most of order 
H. For our purposes the neutrinos can be either Dirac or 
Majorana particles: the details will not affect our consid¬ 
erations. The only necessary ingredient is that, according 
to , the neutrino mass is a function of the scalar held, 
i.e. rrii, = Here we consider three species of neu¬ 

trinos with the same mass and choose a held-dependence 
of the form = Mq exp{/3/), with (3 = 0(1). The form 
of the coupling chosen is well motivated (see e.g. |2,@) 
and has been considered in the past in models with dark 
matter/dark energy interaction [1,I3|- We point out, how¬ 
ever, that results for other potentials and couplings are 

similar to the ones presented here [13- 
In the cosmological context, neutrinos cannot be de¬ 
scribed as a huid. Instead, we must solve the distribu¬ 
tion function /(x/p/r) in phase space (where r is the 
conformal time). We are interested in times when neu¬ 
trinos are collisionless, and so the distribution function 
/ does not depend explicitly on time. Solving the Boltz¬ 
mann equation, we can then calculate the energy density 
stored in neutrinos (/o is the background neutrino distri¬ 
bution function): 

Pi/ = ^ y q‘^dqdnefo{q), (1) 

with = q^+m,^{(j))'^a‘^, a is the scale factor and g* = ap® 
is the comoving momentum. The pressure is 

Pr = ^ J q’^dqdUlfoiq)^. ( 2 ) 


2 


From these equations one can easily derive that 

Pu + ZH{py +pu) = - 'iPi') (3) 

(the dot representing the derivative with respect to r). 
This equation is akin to the equation for matter coupled 
to a scalar field in scalar-tensor theories. The equation 
of motion for the scalar field reads 

ip + 2H^ + a^— = -a^ {p^-3p^), (4) 

which can be obtained from the energy conservation 
equation of the combined fluid of neutrinos and dark en¬ 
ergy 

Pi/ + + ^H(,Pv + P0 + Py + P(j>) = 0 (5) 

and eq. 0 . From this equation and our choice of 
the dynamic of the field is specified by the effective po¬ 
tential 

Kff = V{ct>) + (p, - 3p,)e^^, (6) 

where p^ = pye~^^ and pi, = are independent 

of (j). With V = the field value at the 

minimum of the effective potential is given by pmin = 
A“^ ln(-\/3AVb/-\/2/3(pi/ — 3pi,)). Restricting to the case 
A > 0, the effective minimum only exists for /3 > 0. 

With our choice of potential and coupling, the setup is 
similar to that studied in m For such a system, a num¬ 
ber of critical points have been identified of which only 
two of these are compatible with an accelerating universe. 
However in our work there is the major difference that the 
dark energy field couples to neutrinos rather than cold 
dark matter (CDM). Furthermore neutrinos have never 
dominated the dynamics of the universe in the past and 
their equation of state is not constant. *** When the neu¬ 
trinos are relativistic, the coupling terms containing the 
trace of the neutrino energy momentum tensor are small, 
but non-zero. In particular, — 3p„)/pi, <C 1. Given 
that (j) is at most of order H in quintessence models, and 
[3 is of order one, it is obvious that the second term on 
the left hand side in eq. © dominates over the coupling 
term. We have confirmed this by numerically solving the 
Boltzmann equation. The source terms are only signifi¬ 
cant when the neutrinos become non-relativistic. *** 

Typically, we find that the system passes through a se¬ 
ries of four stages. Firstly, when the neutrinos are ultra- 
relativistic, the field is frozen and the neutrino mass is 
constant. Then, as the neutrinos start to become non- 
relativistic, part of the energy of the neutrinos is trans¬ 
ferred to kinetic energy for the scalar field. Subsequently, 
at a temperature close to the neutrino mass, the neutri¬ 
nos become non-relativistic and begin to scale similarly 
to dark energy but differently compared to dark matter. 
Here the kinetic energy dominates the dynamics of the 
scalar field. During this time the neutrino mass starts 
to evolve significantly. Finally, typically at a redshift of 
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FIG. 1: Background evolution: In the upper panel, we plot 
the evolution of the density parameters for a model with /3 = 
0, A = 1. In the lower panel the corresponding plot with 13 = 1 
is shown. (Neutrinos: solid line, CDM: dot-dashed line, scalar 
field: dotted line and radiation: dashed line.) In all cases, 
the mass of the neutrinos is mi, = 0.314 eV today. We are 
considering a flat universe with = 0.022, Q.ch^ = 0.12, 

Q.i,h^ = 0.01 and h = 0.7. 


order unity, the energy density of dark energy takes over 
and starts to dominate, while the other energy densities 
decay away. 

In Figures 1 and 2 we plot the evolution of the density 
parameter fli = pi/pcr (pcr is the critical density) of the 
matter species for different choices of /3 and A. Clearly, 
the interaction between the scalar field and the neutrinos 
significantly modifies the evolution of the neutrino den¬ 
sity. The differences in the evolution of Ocdm and for 
different choices of (3 and A are obvious. The neutrinos 
do not scale like a~^ for some time, since the pressure 
does not vanish completely at the beginning of the mat¬ 
ter dominated epoch and also because of the coupling 
to the scalar field. This has important consequences for 
the early time integrated Sachs-Wolfe (ISW) effect, as 
we will discuss shortly. 

For the model studied here, the neutrinos are heav¬ 
ier in the past and become lighter during the cosmolog¬ 
ical evolution. This results in a larger neutrino density 
in the past, for the times when the neutrinos are non- 
relativistic. Since the mass of the neutrinos is consider¬ 
ably larger in some models, this affects the freestreaming 
length and hence the matter power spectrum (see Fig¬ 
ure 3). In Figure 2 (lower panel) we have also plotted 
the evolution of the neutrino mass for different choices 
for [3 and A. In models with positive (3 the effective 
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FIG. 2: The upper plot is the same as Figure 1, but choosing 
(3 = —0.79 and A = 1. The cosmological parameters are 
chosen as in Figure 1. The lower plot shows the evolution of 
the neutrino mass in the different models (solid line: /3 = 0, 
A = 1; short dashed line: /3 = 1, A = 1; dotted line: /3 = 
—0.79, A = 1; long dashed line /3 = 1, A = 0.5.) 

potential possesses a minimum, which explains the late 
time increase of the masses for the cases with /3 > 0: 
the field rolls from large field values towards the mini¬ 
mum, overshoots, comes to a halt and is currently rolling 
back towards the minimum. In this case, starting from 
about redshift unity, the neutrino energy density will de¬ 
cay more slowly than in the cases where the minimum 
does not exist. This particular behaviour plays a role in 
the late time ISW effect (see below). 

In order to study the evolution of cosmological pertur¬ 
bations, we have extended the Boltzmann treatment of 
M to include the coupling between neutrinos and dark 
energy ^ (see e.g. 0 for detailed discussions on the 
physics of CMB anisotropies). To calculate the power 
spectra, we modified the CAMB code M accordingly. 
In Figure 3 we plot the anisotropy spectrum for different 
choices of (3 and A. We observe a number of differences 
with respect to the uncoupled case. Firstly, we see an in¬ 
crease in power on scales larger than a degree (multipole 
number I < 100). In some cases an interesting reduc¬ 
tion of power can be observed on larger scales (multipole 
number I < 10). Furthermore, for some choices of pa¬ 
rameters, the positions as well as the relative heights of 
the peaks are also affected. We will now discuss these 
effects in more detail. 

As discussed above, the background evolution is mod- 
ihed in the presence of mass-varying neutrinos. In par¬ 
ticular the density of neutrinos is larger at early times 


FIG. 3: Upper panel: the CMB anisotropy spectrum (nnnor- 
malized). Solid line: /? = 0, A = 1; short-dashed line: 13 = 1, 
A = 1; dotted line: (3 = —0.79, A = 1; long-dashed line: 
/3 = 1, A = 0.5. The lower panel shows the matter power 
spectrum. From the top curve to the bottom curve: {/3 = 0, 
A = 1), (/3 = 1, A = 0.5), (/3 = -0.79, A = 1). The matter 
power spectrum for (/3 = 1, A = 1) is indistinguishable from 
the (/3 = 0, A = 1) curve. 
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FIG. 4: Evolution of the sum of the metric perturbations 
"F -b 5*. Solid line: /3 = 0, A = 1; short-dashed line: 13 = 1, 
A = 1; dotted line: (3 = —0.79, A = 1; long-dashed line: 
(3 = 1, X = 0.5. The scale is k = 10“®Mpc“^. 

in models with P 0. Around the period of matter- 
radiation equality, the coupling of the neutrinos to the 
scalar field causes the neutrino density to decay faster 
than the energy density of CDM (even if the pressure of 
the neutrinos is negligible). This can be seen from eqn. 
© and using the fact that during this period P<p < 0 
for the models under consideration. As a result, the 
regime between the radiation and matter dominated era 
is prolonged, which can be seen in Figures 1 and 2. This 
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implies that the evolution of the metric perturbations 
<^500 = —2a^'I' and 5gij = —2a^^5ij in this period is sig¬ 
nificantly modified, as shown in Figure 4. The integrated 
Sachs-Wolfe effect (ISW) is an integral of 4'4-4) over con¬ 
formal time and wavenumber k and therefore depends on 
the parameters [3 and A. The changes in the evolution of 
4) and $ in the redshift range z « 50 — 1000 imply an 
excess of power in the CMB spectra, which can be seen in 
Figure 3 in the multipole range 10 < / < 100 for models 
with ^ 0. 

The anisotropies on very large scales {I < 20) are domi¬ 
nated by the late time ISW, i.e. by the evolution of T-l-d) 
in the redshift range between z = 0 and z k, 1, which is 
governed by the evolution of the background and the per¬ 
turbations. In particular, and pi, as well as the equa¬ 
tion of state of dark energy affect the late time behavior 
of cosmological perturbations. As mentioned above, the 
evolution of the scalar field is influenced by the presence 
of a coupling to the neutrinos and hence the equation of 
state of dark energy depends upon /3. Likewise, the clus¬ 
tering properties of dark energy depends on the coupling 
to neutrinos (see for a discussion on the clustering 
of dark energy and its impact on the CMB). The neu¬ 
trinos will generally tend to fall into the potential wells 
of dark matter, although at a rate slightly dependent on 
the coupling to the scalar field. The scalar field itself 
will cluster together with the neutrinos and thereby af¬ 
fecting the gravitational potential. For some choices of 
f3 and A we find a suppression of power relative to the 
case with /3 = 0. In particular, the anisotropy spectrum 
for the model with (3 = 1 and A = 1 only differs from 
the uncoupled case on large angular scales, leaving the 
acoustic peaks almost unmodified, and should therefore 
result in an improved fit with the latest WMAP data 
[T^. However, the reduction of power on large angular 
scales is not generic and other choices for (3 and A lead 
to an enhancement of power in the region 1 = 2 — 100, as 
can be seen in Figure 3. 

Finally, the shifts and slight rescaling in the peaks is 
caused by the different densities stored in massive neutri¬ 
nos, baryons and CDM at the time of decoupling when 
changing the parameters (3 and A. The physics is very 
similar to the cases studied in na. The predicted mat¬ 
ter power spectra look very similar to standard models 


with CDM 4- hot dark matter. The damping observed 
in the spectra can be simulated by an averaged neutrino 
mass in the models considered here. However, there are 
new signatures in the CMB power spectra which cannot 
be obtained with an averaged neutrino mass and are due 
to the coupling between dark energy and neutrinos. 

We would like to point out that the decay of neutri¬ 
nos into (^—quanta does not play a role for the parame¬ 
ters chosen here. Potenti ally, this can have an important 
effect in cosmology (see [l^)- The Lagrangian for the 
neutrinos is 

JO-i. = rriuifOvv « Mqvv + ^^(0 - 4- ..., (7) 

Mpi 

where we used m,y{<f>) = Mq exp{(3(j)), expanded around 
the minimum ^min and have neglected higher order 
terms. This Lagrangian has the same form as the one 
used in 0 if we identify the coupling to be = 
f3Mo/Mpi. For (3 = 0{1) and Mq « C>(eV) the cou¬ 
pling 5 <C 1 and indeed much smaller than the value 
[g = 10-5) used in [3. 

In conclusion, cosmologies with neutrino-dark energy 
coupling have a rich phenomenology. It is clear from the 
results presented in this letter that models with mass- 
varying neutrinos cannot be mimicked with an averaged 
constant neutrino mass. We have found that some mod¬ 
els with a coupling of the order of the inverse of the 
Planck mass present a reduction of power in the temper¬ 
ature CMB anisotropies spectrum at low multipoles but 
a standard cosmology peak structure in line with current 
CMB data. Our work implies that CMB anisotropies as 
well as large scale structures will be able to constrain 
parameters of the theory tightly. In future we will inves¬ 
tigate other potentials and couplings, elaborate on the 
degeneracies between the parameters and will use cur¬ 
rent data to constrain such models [Tn|. 

We are grateful to O. Elgaroy, D. Hooper, A. Mel- 
chiorri, J. Silk and C. Skordis for useful discussions 
and to A. Lewis for allowing us to use the CAMB 
quintessence module. The authors were supported by 
PPARC (AWB) and the Research Council of Norway 
through project number 159637/V30 (DFM). DTV ac¬ 
knowledges a Scatcherd Scholarship. 


[1] A.G. Riess et al., Astrophys.Journ. 607 , 665 (2004) 

[2] C. Wetterich, Nucl. Phys. B302, 645 (1988) 

[3] P.J.E. Peebles and B. Ratra, Astrophys. J. Lett. 325, 17 
(1988) 

[4] M. Doran and J. Jacket, Phys.Rev.D 66, 043519 (2002) 

[5] S.M. Carroll, Phys.Rev.Lett. 81 , 3067 (1998) 

[6] C. Wetterich, Astron.Astrophys. 301 , 321 (1995); see 
also G.W. Anderson and S.M. Carroll, astro-ph/9711288 

[7] L. Amendola, Phys.Rev.D62, 043511 (2000) 

[8] R. Pardon, A. E. Nelson and N. Weiner, JCAP 0410:005 
(2004) 


[9] X.-J. Bi, P. Gu, X. Wang, X. Zhang, Phys.Rev.D69, 
113007 (2004); P. Hnng, H. Pas, astro-ph/0311131 
D.B. Kaplan, A.E. Nelson and N. Weiner, Phys.Rev.Lett. 
93 , 091801 (2004); R.D. Peccei, Phys.Rev.D 71 , 
023527 (2005); E.I.Guendelman and A.B. Kaganovich, 
hep-th/0411188 X. Bi, B. Feng, H. Li and X. Zhang, 
hep-ph/0412002 V. Barger, P. Hnber and D. Marfatia, 
hep-ph/0502196 M. Cirelli, M.C. Gonzalez-Garcia and 

C. Pena-Garay, hep-ph/0503028 

[10] A.W. Brookfield, C. van de Bruck, D. F. Mota and 

D. Tocchini-Valentin!; in preparation 





5 


[11] D. Tocchini-Valentini and L. Amendola, Phys.Rev.D 65, 
063508 (2002) 

[12] C.-P. Ma and E. Bertschinger, Astrophys. Journ. 455, 7 
(1995) 

[13] S. Dodelson, Modern Cosmology, Academic Press (2003), 
R. Durrer, J. Phys. Stud. 5, 177 (2001) 

[14] A. Lewis, A. Challinor and A. Lasenby, Astrophys. Journ. 
538, 473 (2000) 

[15] J. Weller and A. Lewis, MNRAS 346, 987 (2003); R. 


Bean and O. Dore, Phys.Rev.D 69, 083503 (2004) 

[16] C. L. Bennett et ah, Astrophys. J. Suppl. 1, 1 (2003) 

[17] S. Dodelson, E. Gates and A. Stebbins, Astrophys. J. 
467, 10 (1996) 

[18] J.F. Beacom, N.F. Bell and S. Dodelson, 

Phys.Rev.Lett.93, 121302 (2004); S. Hannestad, 

JCAP0502:011 (2005) 



